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ABSTRACT 

Modifications of hadron masses and some of their basic properties with 
temperature or nuclear density are considered as one possible signatures for the 
formation of dense hadronic matter in nuclear collisions. We discuss here some 
basic results on the temperature dependence of hadron properties obtained 
from calculations of hadron correlation functions, the chiral condensate as well 
as the equation of state in finite temperature lattice QCD. 



1. Introduction 

The theory of strongly interacting particles - Quantumchromodynamics (QCD) 
- has many unusual non-perturbative features, which are under intensive theoret- 
ical and experimental investigation. It is the hope that some of the low energy 
characteristics of QCD like confinement and chiral symmetry breaking may become 
better understood through studies of strongly interacting matter, ie. large sys- 
tems of hadrons at high temperature and/or high baryon number density. This 
is, of course, also of relevance for our understanding of the evolution of the early 
universe. 

Properties of strongly interacting matter are studied in ultra-relativistic heavy 
ion experiments. Through an analysis of the final state hadrons one tries to extract 
information about the dense thermal medium, created in these collisions. Basic 
properties of hadrons like their masses, widths and decay constants reflect many 
features of the complex, non-perturbative structure of the QCD vacuum. Any mod- 
ifications of these hadronic properties at high temperature and/or density may lead 
to modifications in the particle spectra, which will lead to experimentally observable 
effects. It thus is important to understand theoretically the equilibrium properties 
of hadrons in dense/hot nuclear media, although the detailed description of heavy 
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ion experiments will still involve further assumptions about non-equilibrium effects 
and the thermal evolution in a heavy ion collision 1 . 

Effective theories, deduced from QCD in order to describe low energy proper- 
ties of hadrons, establish a close link between hadronic properties and the non- 
perturbative structure of the QCD vacuum, which is described by various non- 
vanishing condensates. This is, for instance, used in the operator product expan- 
sion (OPE) for correlation functions of hadronic currents 2 , which allows to relate 
hadron masses to condensates of the quark and gluon fields. The temperature and 
density dependence of the latter has been discussed within the context of chiral 
perturbation theory 3 . Through an extension of the OPE to finite temperature it is 
then possible to discuss the temperature dependence of hadron masses and other 
hadronic properties in a nuclear medium. We should, however, mention that this 
approach is not at all straightforward and led to some controversial discussion about 
the relevant states which should be used in the finite temperature expansion 4 . 

We will concentrate here on a discussion of the influence of non-zero tempera- 
tures on hadronic properties^. The applicability of the OPE as well as chiral per- 
turbation theory is limited to low temperatures. In the vicinity of the QCD phase 
transition from hadronic matter to the quark gluon plasma phase, however, the 
hadronic medium becomes so dense that the influence of many particle states and 
resonances can no longer be neglected. In this regime an entirely non-perturbative 
approach becomes mandatory. This is the regime where lattice QCD might provide 
the only reliable approach to study the temperature dependence of condensates, 
hadrons and their properties. 

We will discuss here some of the recent results of lattice calculations, which aim 
at an analysis of hadronic properties close to the QCD phase transition. In the next 
section we briefly review results about the thermodynamics of QCD, focusing on 
the equation of state and chiral properties in the vicinity of T c . Section 3 is devoted 
to a discussion of lattice calculations of hadron masses and decay constants at finite 
temperature within the quenched approximation. We briefly comment about the 
situation in two-flavour QCD in Section 4. Finally we give our conclusions in Section 
5. 

2. Basic Thermodynamics from Lattice QCD 

Thanks to the rapid development of computer technology and the improvement 
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of numerical algorithms for Monte Carlo calculations the studies of lattice QCD at 
finite temperature lead to a steadily improving understanding of details of the QCD 
phase transition. In the absence of dynamical quark degrees of freedom (pure gauge 
theory/quenched QCD) we have obtained quantitative results for the order of the 
transition, the transition temperature and equation of state on lattices of different 
size. This allows with some confidence an extrapolation to the continuum limit 5 ' 6 . 
Although simulations with light quarks did not yet reach a similar accuracy we 
know also here a lot about the qualitative behaviour of QCD at finite temperature. 

As discussed in the introduction we expect that the properties of hadrons at 
finite temperature are closely related to changes in non-perturbative condensates in 
the QCD vacuum. With increasing temperature more and more hadronic states get 
excited and populate the vacuum. This is expected to result in decreasing values for 
the condensates. Finally, when the density of hadrons becomes too large the whole 
space is filled by hadronic bubbles which completely suppress the non-perturbative 
chiral condensate - the phase transition to the quark gluon plasma phase occurs 
and the chiral symmetry is restored above T c . 

This simplified picture of the QCD phase transition suggests a similarity with 
a percolation phase transition. In fact, the close relation between the density of 
the hadronic medium, a percolation threshold, and the occurrence of the QCD 
phase transition finds support in the numerical simulations of QCD with varying 
number of partonic degrees of freedom. With increasing number of light (nearly 
massless) quark flavours the phase transition has been found to occur at lower 
temperatures. This seems to be natural for a phase transition controlled be a 
percolation threshold: QCD with n/ massless quarks describes a world with (rii — 1) 
Goldstone particles (pions). At a given temperature the density of the hadronic 
medium thus will be proportional to (n* — 1) and a critical density at which these 
hadrons start overlapping occurs earlier with larger rif. While the change in the 
critical temperature is relatively small when changing the number of flavours from, 
in general, two to four, there is a large change between n/ = (pure gauge theory) 
and nf = 2. The transition temperature in QCD with two light quarks is found 
to occur at T c (rif = 2) ~ 150MeV, while in the purely gluonic theory, nj = 0, it 
is found to be substantially higher, T c (rif = 0) ~ 230MeV. This may also easily 
be understood in terms of a percolation threshold: In the absence of quarks there 
are no light hadrons, the lightest states are glueballs with a mass of 0(1 GeV). 
One thus needs a rather large temperature to excite enough of these heavy glueball 
states to reach a critical density. 



Besides this change in the temperature scale and details of the phase transition 
itself - the transition is first order for rif = and likely to be second order for 
nf = 2 - the temperature dependence of bulk thermodynamic quantities is very 
similar. In Fig.l we show recent results for the energy density, e, and the pressure, 
P, in the SU(3) gauge theory 5 (n/ = 0) as well as two-flavour QCD 7 . As can 
be seen, in both cases a substantial change in the effective number of degrees of 
freedom occurs only in the vicinity of T c . 



Figure 1. Energy density and pressure for the pure SU(3) gauge theory (a) and two-flavour 
QCD (b) as a function of T/T c . The pure gauge theory results have been obtained from simulations 
on a large 32 X 6 lattice while those for two-flavour QCD so far exist only for a rather small 
lattice (8 X4). For more details on the inherent finite size effects in these lattice calculations we 
refer to the above papers as well as one of the reviews on finite temperature lattice QCD . 

The numerical results for the energy density suggest that up to temperatures of 
about 0.9T C the hadronic medium is well described by a gas of hadrons, in which 
only the lightest states contribute. In the case of a SU(2) gauge theory the critical 
energy density at T c has been studied in great detail 8 . Here it is found that the 
deconfinement transition occurs at e c /T^ ~ 0.25, which is only a factor five larger 
than the energy density of a glueball gas. Already below 0.9T C higher excited states 
play no role and one is left with a rather dilute glueball. A rough estimate yields 
values of about (0.1-0.2) for the number of glueballs per hadronic volume. 



The slow variation of the energy density and pressure below 0.9T C shows that also 
the gluon condensate will vary only little with temperature. The gluon condensate 
is directly related to the difference e — 3P, 

(G 2 ) T = (G 2 ) -(e-3P) , (1) 

where (G 2 )q denotes the gluon condensate at T = 0. The small variation with T, 
which can be deduced from Fig.l, is in accordance with chiral perturbation theory, 
which suggests that the temperature dependence of (G 2 )t starts only at 0(T 8 ) 
(Ref. 3). 

For the chiral condensate (ifjifj) chiral perturbation theory suggests a stronger 
dependence on temperature, A 

T n f - 1 ( T 2 \ n 2 f - 1 / T 2 \ 2 J 



Figure 2. The chiral condensate at finite temperature extrapolated to the chiral limit 
(rriq =0) and normalized to the corresponding zero temperature value for various values of the 
number of flavours (Fig. 2a). In Fig. 2b we show results from a simulation in quenched QCD (nj=0) 
for various values of the quark mass. The simulations have been performed at fixed value of the 
gauge coupling , 6/g =6.0, on a 32 X8 lattice, corresponding to T~0.92T C and on various low 
temperature lattices with sizes ranging from 24 X32 to 32 X64 (Ref.10). The condensates (lplp)i 
and \1p1p)2 are defined in Section 4. No temperature dependence is visible. 



In Fig. 2 we show a collection of results obtained for QCD with various num- 
ber of flavours. They suggest at most a weak dependence of the chiral condensate 
on temperature. However, in particular in the case of QCD simulations with light 
quarks this behaviour has to be confirmed in simulations on larger lattices. In the 
case of pure SU(3) gauge theory simulations on large lattices have been performed 9 
and confirm that there is no significant temperature dependence of the chiral con- 
densate up to T ~ 0.92T C . It thus seems to be conceivable that significant changes 
in the QCD condensates occur only quite close to T c . 

3. Hadronic Properties close to T c 

Information about hadron masses and decay constants can be extracted in lattice 
simulations from the long-distance behaviour of correlation functions of hadron 
operators 

G H (x) = (H(x)H+(0)) ^ e - mH ^ , x=(t,x) , (3) 

with H(x) denoting an operator with the appropriate quantum numbers of the 
hadronic state under consideration, for instance H (x) = ^ui^l^dix) f° r the pion. 
At zero temperature one studies the behaviour of the correlation function at large 
Euclidean times r. From the exponential decay of Gh one deduces the hadron 
mass of the lightest state in this channel, whereas the amplitude is related to the 
decay constant of this hadronic state. At finite temperature the Euclidean extent 
is limited, < r < 1/T, and one thus studies the behaviour of Gh for large 
spatial separations, \x\ — ► oo. This yields information about the finite temperature 
screening masses which are related to pole masses as long as there is a bound state 
in the quantum number channel under consideration. 

Finite temperature screening masses have been studied in lattice simulations 
of QCD for quite some time 11 . A very drastic qualitative change in the screening 
masses is seen when one crosses the QCD transition temperature. Parity part- 
ners become degenerate above T c , the pseudo-scalar mass becomes massive and 
approaches the free quark/anti-quark value, m meS on = 27rT, at large temperatures. 
These features do not seem to depend much on the number of quark flavours. In 
particular, they also have been found in quenched QCD simulations. It thus seems 
to be meaningful to first study the behaviour of hadronic properties in the quenched 
approximation where results on large lattices with high accuracy can be obtained. 
In the following we will describe the results of such an investigation performed on 



a rather large lattice (32 3 x 8) close to the phase transition temperature. 

3.1. The GMOR Relation and the Pion Decay Constant 

The GMOR relation relates the chiral condensate to the pion mass and pion 
decay constant, 

flml = m q ($ijj) mq =o ■ (4) 

Below T c the pion is a Goldstone particle, its mass squared depends linearly on the 
quark mass, 

ml = a^rriq . (5) 

A calculation of the chiral condensate and the pion mass at different values of the 
quark mass allows the determination of the pion slope, a^, and the zero quark mass 
limit of the condensate, (ipip) m =o. 

The pion decay constant f n can be determined directly from the relevant matrix 
element 

V2f n ml = (01^75^1^+) . (6) 

The square of the matrix element appearing on the right hand side of Eq.(6) 
is proportional to the amplitude of the pion correlation function and can be deter- 
mined in a Monte Carlo calculation. A comparison of f n determined this way with 
the ratio {^^)m q -o/a^ thus provides a direct test of the GMOR relation at finite 
temperature. This is shown in Fig. 3a at T ~ 0.92T C . Moreover, we can compare the 
value of f n calculated at finite temperature with corresponding zero temperature 
results. This is shown in Fig. 3b. Clearly we do not have any evidence for violations 
of the GMOR relation nor for a significant change of f n with temperature below 
T c . The sudden change above T c reflects the drastic change in the structure of the 
pseudo-scalar correlation function. It does no longer have a pole corresponding to 
a Goldstone-particle (pion) . The existence of a pseudo-scalar bound state above T c 
thus is questionable. Certainly for large temperatures such a state does not exist, 
the correlation function is dominated by a quark/anti-quark cut. The notion of f n 
used for the square root of the amplitude of the correlation function should thus be 
used with caution above T c . 

3.2. Vector Meson Mass and Nucleon Mass 

The possibility of a variation of the p meson mass with temperature has been 
discussed a lot as this might lead to modifications of dilepton spectra 1 , which are 
experimentally detectable. The temperature dependence of the meson masses has 



Figure 3. Test of the GMOR relation at T=0.92T C in quenched QCD (a) and the temperature 
dependence of the pion decay constant (b). In Fig. 3a we compare the results obtained from the 
amplitude of the pion correlation function through Eq.(6) (circles) and extrapolated to zero quark 
mass with the result obtained from the GMOR relation (star). Also shown is the zero temperature 
result at this value of the gauge coupling (square). 

been discussed within the framework of the OPE. Arguments have been given that 
the meson masses are temperature independent up to 0(T 2 ) (Ref. 4). The Monte 
Carlo calculations of the vector meson correlation function at finite temperature 
also show no significant temperature dependence of the mass even close to T c . In 
Fig. 4 we show the result of our calculation of the screening mass in the vector 
channel correlation function at T ~ 0.92T C and compare this with zero temperature 
calculations at the same value of the gauge coupling. There is no evidence for 
any temperature dependence. The same holds true for the nucleon mass, although 
the details are more subtle in this case. As can be seen in Fig. 4b there is a clear 
difference between the local masses, rriN(z) ~ 1uGn(z)/Gn(z + 1), extracted on a 
large zero temperature lattice and those extracted at finite temperature on a lattice 
of size 32 3 x 8. However, as the nucleon is a fermion, there is a non-negligible 
contribution from the non-zero Matsubara mode to the nucleon screening mass, 

m N = ^m 2 N + (ttT) 2 (7) 

After removing the thermal contribution, ttT, the nucleon mass, m^, agrees with 



Figure 4. The p-meson mass in quenched QCD for various values of the quark mass (a) and 
local masses from the nucleon correlation function (b). The simulations have been performed at 
fixed value of the gauge coupling , 6/g =6.0, on a 32 X8 lattice corresponding to T^i0.92T c and 
on various low temperature lattices (see Fig. 2). In Fig. 4b we show estimates for the nucleon mass 
obtained from the nucleon correlation function at distance z. These local masses converge to the 
nucleon mass in the limit z — ► OO. The horizontal band indicates the corresponding extrapolation 
at zero temperature. No temperature dependence is visible in the p-meson mass as well as the 
nucleon mass (see text for appropriate subtraction of the finite temperature Matsubara mode). 

the zero temperature result within statistical errors. 

4. Two-Flavour QCD 

The discussion of hadron masses and decay constants in quenched QCD shows 
that there are no visible deviations from the zero temperature values up to tem- 
peratures T ~ 0.9T C . One may expect that this is partly due to the fact that the 
phase transition in pure £77(3) gauge theory is first order. The change from the 
hadron phase to the quark gluon plasma phase is more abrupt than in the case of 
a second order transition and might lead to less significant changes in condensates 
and hadronic properties. 

The presently available simulations of QCD with two light quark flavours, in- 
deed, suggest that the transition is smoother. In fact, the analysis of the scaling 
behaviour of various thermodynamic quantities with the quark mass seems to be 



consistent with the behaviour expected for a second order chiral phase transition 12 . 
We thus may wonder whether stronger modifications of hadron properties can be 
expected in this case. 



Figure 5. The chiral condensate (a) and the chiral susceptibility (b) in two flavour QCD 
obtained from simulations on an 8 X4 lattice 

In Fig. 5 we show the behaviour of the chiral condensate and its derivative with 
respect to the quark mass - the chiral susceptibility % m , 

Xm= z^ z ' (8) 

V <9m 2 

Also shown there is an improved estimator for the zero quark mass chiral condensate 
in which the 0(m q ) contribution to (V'V')i * s removed, 

(4>ip)2 = (4>ip)i - m q Xm ■ (9) 

As can be seen (V>V>)i deviates from ('0V') 2 strongly only in a narrow region around 
T c . This indicates that only in this region contributions from the singular part of 
the free energy, which is responsible for the occurrence of a phase transition, is 



dominant. Also the chiral susceptibility, % m , has a very narrow peak at the pseudo- 
critical point. We thus expect that also in the case of two-flavour QCD substantial 
changes in hadronic properties will occur only in a narrow temperature regime close 

tO T c (m q ) = Tpeak, 



^crit 



J- peak/2 -'peak 



■ peak 



<0.1 . (fO) 



5. Conclusions 

We have discussed some of the basic properties of QCD thermodynamics, which 
are relevant for a discussion of the temperature dependence of hadronic properties 
in the hadronic phase of QCD. There are strong indications that the gluon as 
well as the chiral condensates show at most a weak dependence on temperature 
for T<0.9T C . So far no statistically significant temperature dependence could be 
observed in lattice QCD calculations. The same is true for hadron masses and decay 
constants calculated in quenched QCD. 

At least for the chiral condensate one expects, however, a strong variation with 
temperature for 0.9T C < T < T c . This is supported by the currently available lattice 
calculation in two-flavour QCD. The temperature variation of the gluon condensate, 
on the other hand, might not be that strong as indicated by the behaviour of the 
equation of state shown in Fig.l. The influence of this on the behaviour of hadronic 
parameters will certainly be investigated in the near future. 
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